A not-so-normal mode decomposition 
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We provide a generalization of the normal mode decomposition for non-symmetric or locality constrained 
situations. This allows for instance to locally decouple a bipartitioned collection of arbitrarily correlated oscil- 
lators up to elementary pairs into which all correlations are condensed. Similarly, it enables us to decouple the 
interaction parts of multi-mode channels into single-mode and pair-interactions where the latter are shown to be 
a clear signature of squeezing between system and environment. In mathematical terms the result is a canonical 
matrix form with respect to real symplectic equivalence transformations. 



The normal mode decomposition is a ubiquitous and indis- 
pensable tool in physics and engineering. It allows us to trans- 
form into frames where seemingly complex systems decouple 
into elementary units each of which can be tackled individu- 
ally. Consider for instance a collection of harmonically cou- 
pled oscillators, classical or quantum. Then the normal mode 
decomposition provides us with a canonical transformation 
which decouples the state of all oscillators into independent 
normal modes. Similarly, by the seminal work of Williamson 
lUl], the evolution of a system governed by a quadratic Hamil- 
tonian can be decoupled into independent elementary parts. In 
both cases — states and evolutions — this decoupling requires 
a global transformation. In many situations, however, addi- 
tional restrictions prevent us from simplifying the structure of 
a system by such a transformation. We might for instance 
want to preserve the amount of correlations or entanglement 
of a bipartitioned collection of oscillators. Or, in the case of 
evolution, we might be faced with dissipation, i.e., an inacces- 
sible environment. 

In this work we provide a generalization of the normal 
mode decomposition which is applicable in these and other 
situations and allows again to decouple correlations and inter- 
actions into elementary parts. This will reveal a remarkable 
structure of both, correlations and interactions, in which ele- 
mentary units turn out to be pairs of modes rather than single 
modes. Hence, as for the normal mode decomposition, we 
get again a significant simplification in many contexts dealing 
with many-body harmonic or bosonic systems, albeit with a 
richer structure. 

In the first part we will state and prove the main result in 
terms of matrix analysis where it amounts to a normal form 
with respect to symplectic equivalence transformations, i.e., 
a symplectic analogue of the singular value decomposition — 
inspired by recent advances in symplectic geometry iSB- In 
the second part we will then apply it to the above mentioned 
cases of bipartite correlations and dissipative evolutions (de- 
picted in Figs. 1,2) and connect it to known results. 



PRELIMINARIES 

Before introducing some basic notions let us mention that 
although we will have quantum systems in mind in the fol- 
lowing, all results hold for classical systems in exactly the 



same way. Similarly, note that the oscillators do not have to 
be mechanical but might as well correspond to electromagnet- 
ical field modes, charge-phase oscillations in Josephson junc- 
tions or collective spin fluctuations, e.g., in atomic ensembles. 
In each of these cases we have pairs of canonically conjugate 
variables obeying the same (commutation) relation as position 
and momentum. 

Consider now n quantum mechanical oscillators char- 
acterized by a set of momentum and position operators 
(Pi, . . . , P„, Qi, . . . , Qn) ='■ R which obey the canonical 
commutation relations [Rk, Ri] = icFku with 
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the symplectic matrix. A canonical/symplectic transformation 
maps Rk 1-^ X)/ SkiRk such that the commutation relations 
(or classical Poisson brackets) are preserved, i.e., SaS'^ = a. 
We will denote the group of real symplectic transformation on 
n modes by Sp{2n). For a basic introduction into symplec- 
tic transformations and the appearance and use of canonically 
conjugate variables in quantum information theory we refer 
the reader to 13] and 101. 

The essence of the ordinary normal mode decomposition 
is the fact that for any positive definite matrix X G M^" ^ ^" 
there is an G Sp{2n) such that 



SXS^ = diag(i/i, 



(2) 



In case X represents a Hamiltonian H = J^ki -^kiRkRi the 
i/k are the normal mode frequencies. If Xki = {{Rk — 
(Rk), Ri — {Ri)}+) is a covariance matrix, then {uk — l)/2 
is the mean occupation number (phonons/photons) in the fc'th 
normal mode. 



CANONICAL FORM 

We aim at deriving a canonical form for general (not nec- 
essarily symmetric) matrices X E ]R2nx2n ^j^jg]- symplectic 
equivalence transformations X i-^ SiXS2- To this end we 
will first construct a set of invariants which play a role similar 
to the normal mode frequencies i/k- 
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Proposition 1 (Invariants) The eigenvalues of := 
XaX^ (7^ are invariant with respect to symplectic equiva- 
lence transformations of the form X i— > S\XS2- 

The proof of this statement is simple. We have just to exploit 
that SicrSf = a and that for any two matrices AB and BA 
have the same non-zero spectrum. 

Note that the entire spectrum of is two-fold degen- 

erate, i.e., speciT.) = {Ai, . . . , A„, Ai, . . . , A„} 0] and in 
addition complex eigenvalues come in conjugate pairs A, A 
(as it holds for every real matrix). If X is positive definite 
we recover, in fact, the normal mode frequencies Vk = \/Afe 
so that all invariants are positive real numbers in this case. 
The occurrence of complex eigenfrequencies is, in fact, a well 
known phenomenon in many fields of physics in particular 
where harmonic approximations are used in intermediate en- 
ergy regimes. Examples can be found in contexts from molec- 
ular condensates to gravitational waves ^ and sonic black 
holes H. 

While Prop [T] evidently holds for arbitrary rectangular and 
possibly singular matrices X, we will for the sake of simplic- 
ity restrict ourselves to non-singular square matrices in the 
following. Our aim is to show that Ai , . . . , A„ are the only in- 
variants and that they essentially determine the normal-form 
of X with respect to symplectic equivalence transformations: 

Proposition 2 (Canonical form) For every nonsingular ma- 
trix X G M2nx2ri f/jg^g exist real symplectic transformations 
Si , 5*2 such that 



S1XS2 







J 



Ji(Ai] 






J2(A2) 



(3) 

where each Jk is a real Jordan block corresponding to ei- 
ther a complex conjugate pair {A, A} of eigenvalues ofY,{X) 
or to one of its real eigenvalues. In the former case the di- 
agonal of Jk{^) is build up out of real 2x2 blocks of the 
form 



-b a 



X ^ a + ib 



(4) 



Proof. First note that is a skew-Hamiltonian matrix, 

i.e., (Scr)"^ = — (Scr) (we drop the dependence on X in 
the following). For every real skew-Hamiltonian matrix there 
exists a real symplectic similarity transformation such that 
ST,S-^ = -{M e M'^) is block diagonal |l2t|. Exploiting 
in addition that every real matrix M can be written as a prod- 
uct of two real symmetric matrices AB = M 11 Oil we can 
write 



iSX)a{SXfa 



A 
B 



(5) 



Now define a transformation S' by imposing SXS' = W. In 
fact, S' is symplectic as can be seen from 

S'^aS' = W^a'^ [<j{SX)-^a{SX)-^] W ^ a, (6) 



where we have used that a = 1, = aW and that 

the expression in squared brackets is by Eq.Q equal to W~^. 

cr^ = aW is easily seen by exploiting the A, B-block 
structure of W. Note that Eq.® is the point in proof where 
we use non-singularity of X. 

To proceed we exploit the subgroup GL{n) C Sp{2n), i.e., 
the fact that for every real invertible matrix G, the block ma- 
trix G^^ ® G'^ is symplectic Multiplying SXS' from 
the left with A~^ and from the right with a, which is 
a symplectic transformation in its own right, we obtain 1„ 
{—A^B). This can be brought to the claimed form in Eq.(|3]l 
via a symplectic similarity transformation by G^^ G. Here 
we use the real Jordan canonical form J = G{—A^B)G~'^ in 
which complex conjugate pairs of eigenvalues correspond to 
real 2x2 matrices of the form in Eq.(|4|i (cf. lfTlll ). It remains to 
show that the spectrum of is a doubling of the spectrum 
of J. By Proplllwe have that spec{Y.{X)) = spec{T.{l J)) 
so that the identity I](l0j) = j'^ (B J completes the proof. 

□ 



Some remarks on the normal form in Eq.© are in order 
First note that it is minimal in the sense that the number of 
continuous parameters cannot be further reduced by symplec- 
tic equivalence transformations as they are all invariants due 
to ProplT] Similarly looking at the invariants tells us that an 
entirely diagonal normal form, analogous to the usual normal 
mode decomposition, cannot exist in general as real diagonal 
matrices have only real invariants. Hence, there is no way 
of diagonalizing the remaining 2x2 blocks since they cor- 
respond to complex A's. Concerning the possible appearance 
of defective parts in the Jordan blocks 19] we note that, as 
usual, they are not stable with respect to perturbations. The 
fact that matrices with non-defective normal form are dense 
can be seen by noting that for any X G ]g2nx2n (jjgj-g exists 
an X^ arbitrary close to it such that spec(S](Xc)) is only two- 
fold degenerate (and non-singular). The corresponding J has 
no degeneracy and is thus non-defective. 

The presented canonical from can be regarded as a gener- 
alization of the seminal results by Williamson jll] and its ex- 
tensions I12II on normal forms of symmetric (not necessarily 
positive) matrices under symplectic transformations. For pos- 
itive definite matrices X the canonical form in Eq.([3l) and the 
usual normal mode decomposition in Eq.© coincide up to a 
simple squeezing transformation and the invariants are related 
via A I, = ly?. 



In the following we will discuss two applications of the 
above result in contexts where the two sets of modes on which 
Si and S2 act on either correspond to two different parties (Al- 
ice and Bob, say) or to input and output of a quantum channel. 
Note that in both cases the two sets need not be of the same 
physical type. A prominent example of that form is a set of 
light modes coupled to modes of collective spin fluctuations 
in atomic ensembles 11 1311 . 
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FIG. 1 : Left: let Alice and Bob each have a collection of oscillators 
with arbitrary correlations between them (indicated by dashed lines). 
Right: there exist local canonical transformations which condense 
the correlations into elementary units (single modes or pairs) and 
eliminate all cross-correlations. 
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FIG. 2: Left: A multi-mode Gaussian channel consists out of an in- 
teraction part X, which directly couples the input modes, and an 
addition of input-independent noise Y . Right: The interaction part is 
decoupled into single-mode and two-mode parts by applying canoni- 
cal transformations before and after the channel. Two-mode interac- 
tions corresponding to a complex invariant A only occur if the global 
system-plus-environment evolution involves squeezing. 



DECOUPLING OF INTERACTIONS 



DECOUPLING AND CONDENSING CORRELATIONS 

Consider a bipartitioned collection of n + n modes as in 
Fig.l. The covariance matrix F e ]]j4rix4)i be parti- 

tioned into blocks 



(7) 



where F^i, Tb are the local covariance matrices and X de- 
scribes correlations between the two parts. A local symplectic 
transformation Fi— > {Sa®Sb)T{Sa®Sb)'^ transforms the 
correlation block as X i-^ SaXS^. Hence, Prop|2]can be 
directly applied to decouple and condense the correlations. In 
the generic case of non-singular, non-defective X we are then 
left with correlated single modes whose correlations are char- 
acterized by real invariants A and with correlated pairs coitc- 
sponding to complex A with a correlation block of the form 
in Eq.dlli. All the cross-covariances between opposite subsets 
which do not correspond vanish (see Fig.l). 

A particular known instance of this result is the case of 
pure Gaussian quantum states. For pure bipartite states cor- 
relations, which are then due to entanglement, and local spec- 
tral properties determine each other This is the content of 
the Schmidt decomposition which in terms of the covariance 
matrix impUes that F^, F^ and X can then be simultane- 
ously diagonalized by local symplectic transformations such 
that Vk ~ \/\ — \k (Afc < 0) is the mean particle number in 
the fc'th normal mode of each site lfl4l[l9ll . 

This pure state normal form simplified investigations in var- 
ious directions like security proofs in quantum cryptography 
1 1511 or the transformation 11611 . localization |17] and charac- 
terization ifisll of entanglement. Prop|2] now provides the 
analogous normal form for mixed states, which again can be 
regarded as a condensation of correlations. As such, it might 
be a useful first step in quantum information protocols which 
use correlations as a resource. 



Let us now consider dissipative evolutions of multi-mode 
quantum systems. An important class of such evolutions are 
those where system plus environment undergo a global canon- 
ical transformation. For the system (with traced out environ- 
ment) this leads to so-called Gaussian or quasi-free channels 
|[l9[ I20I1 realized by optical fibres or, if the transformation is 
in time rather than in space, by quantum memories build upon 
atomic ensembles 12111 . 

Gaussian channels are characterized by a pair of matrices 
X,Y £ r2"x2« satisfying the constraint iX^aX + Y >ia. 
The covariance matrix evolves then according to 



F 



X'^TX 



Y. 



(8) 



That is, X can be regarded as characterizing direct interac- 
tions between the modes and F is a noise-term which is input- 
independent. If we now apply a symplectic transformation be- 
fore and after the channel then X,Y ^ S1XS2, S'jFS'2. We 
can thus again exploit Prop|2]in order to simplify the structure 
of the interaction part of the evolution. This way of encod- 
ing and decoding information sent through the channel has 
been successfully exploited in the context of channel capaci- 
ties of single-mode channels 12211 for which it leads to a sim- 
ple normal form 12311 . For multi-mode channels it allows us, 
in the generic case, to reduce X to two-mode interactions of 
the form in Eq.(|4l) and single-mode parts (see Fig. 2). 

The appearance of pair interactions corresponding to com- 
plex invariants A is, in fact, a signature of a non-number pre- 
serving system-environment interaction. In order to see this 
recall that a global number-preserving transformation can be 
written as 



S 



C D 

-D C 



(9) 



where C + iZ? is a unitary and the block structure in Eq.© 
refers to a decomposition of phase space into position and mo- 
mentum space (rather than system and environment). For the 
reduced system evolution this leads to an X which has the 
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same structure as S but without the restriction of the matrices 
being real and imaginary parts of a unitary. Let us denote the 
corresponding blocks in X by c and d and calculate 



F + cc^ dc^ - 
cd^-dc^ dd^ 



(10) 



As is Hermitian it has indeed only real eigenvalues Afe 

which shows that pair interactions in the normal form which 
correspond to complex A's witness a squeezing-type interac- 
tion between system and environment. 
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